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ABSTRACT 

This volume contains reports of NASA sponsored studies 
in the area of space flight and guidance theory implementation. 
The studies are carried on by several industrial companies. 

This report covers the period from initiation of the studies 
until September 30, 1963. The technical supervisor of the 
contracts is W. E. Miner, Deputy Chief of the Astrodynamics 
and Guidance Theory Division, Aero-Astrodynamics Laboratory, 
George C. Marshall Space Plight Center. 
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SUMMARY 

This volume contains reports of NASA sponsored studies, 
in the area of space flight and guidance theory implementation. 
The studies are carried on by several industrial c om panies. 

report covers the period from initiation of the studies 
until September 30, 1963. The technical supervisor of the 
contracts is W. E. Miner, Deputy Chief of the Astrodynamics 
and Guidance Theory Division of the Aero-Astrodynamics 
Laboratory, George C. Marshall Space Flight Center. 


INTRODUCTION 

This is the first of a series of reports dealing with, 
the implementation of theory being developed by . Astrodynamics 
and Guidance Theory Division of Aero-Astrodynamics Laboratory 
and the associated contractors. 

The term Progress Report No. 1 (2, 3 » will be used 
for "Progress Reports on Studies in the Fields of Space 
Flight and Guidance Theory. The term Implementation Report 
will be used for "Implementation Repor t No . 1 on Studies. in 
the Fields of Space Flight ana guidance Theory . " These terms 
Will be used for reference to the two companion series of 
reports . 

Two problems are presented in this report. The first 
is that of determining feasibility of adaptive guidance for 
lunar orbital rendezvous ("LOR") type missions. The second 
is that of determining performance data as a basis for . 
evaluating adaptive guidance in non-catastrophic abort. from 
LOR type missions. Both works are preliminary and limited 
in scope. 
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Two papers are presented on the study of feasibility 
of adaptive guidance for LOR type missions. The first is 
by R. S. Polovitch and W. B. Morgan of Boeing Company. The 
second paper is by Dr. S. Hu of Northrop Corporation. The 
problem was arbitrarily defined by MSFC. No attempt was 
made to include the performance problem. Rather, the 
definition was made so that small maneuvers in yaw were 
required. This checked the guidance capability. Boeing 
was assigned the problem of launching into a space-fixed 
conic. This conic was defined as having circular velocity 
of a specified altitude with a 90° path angle (measured 
from vertical). It is realized that due to the earth's 
obliqueness, this will not define a specified conic. Prom 
the practical point of view, the simplified definition of 
the orbit was sufficient. Prom the point of injection into 
the circular orbit, Northrop Corporation developed the 
guidance to insertion into the lunar transit. This latter 
orbit was defined using the JPL lunar deck. Because of 
certain assumptions therein, the resulting lunar transits 
do not make up a continuous group. This caused large errors 
at the moon while having small errors compared to the defined 
end-conditions . 

It is planned in the next implementation report to 
rework the problem utilizing a newly developed earth-moon 
deck and also to elaborate on performance considerations. 

In addition, various engineering constraints will be considered. 
It is also planned to extend the launch window study with a 
view toward determining whether or not additional holds on 
the launch pad are possible. Improvements in procedures will 
be incorporated, especially in the area of curve-fit techniques. 
Progress Reports 3 and 4 refer to articles in this area. In 
these two reports, many ideas are being checked and "packaged" 
for computer use. Lastly, experience alone will greatly 
improve the next results. 

Two papers are presented on the study of performance data 
as a basis for evaluating adaptive guidance for non-catas trophic 
abort from LOR type missions. The first paper is by V. V. Moore 
and P. G. Bourque of Boeing Company. The second paper on reentry 
is by Ann Muzyka and H. Elmore Blanton of Raytheon Company. 

These two studies are complementary. The Raytheon study should 
define acceptable reentry conditions and maximum and minimum 
ranges for free flight through the atmosphere. These reentry 
conditions then become the desired end-conditions for the 
Boeing work. The parameters to be matched here are path angle 

velocity (v), and altitude (y or h). The last parameter 
(h) has been arbitrarily frozen and will only be opened at 
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a very late date. Several points of mismatch can be seen 
in the work. These will be corrected in future work. It 
may be noted that even for preliminary evaluation both papers 
should be considered simultaneously. Both papers offer modi- 
fications and extensions. Many of these will be followed 
for future work. Not clearly fetated is the fact that in 
many cases return to some specified orbit and a later return 
to the earth is possible. This will greatly reduce the size 
of the required landing area. Future work will expand as 
discussed above and then add the full guidance problem. 

Future implementation reports will add the efforts of 
several other contractors. 
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PATH-ADAPTIVE GUIDANCE FOR SATURN V 
THREE-DIMENSIONAL ASCENT TO ORBIT 






by 

R. S. Polovitch 
W. B. Morgan 


SUMMARY 


A 


Pitch and yaw guidance polynomials for the Saturn V vehicle 
have been generated. These polynomials are designed to steer the 
vehicle into a space-fixed orbit at an altitude of 100 n.mi. The 
polynomials are capable of acceptable orbit injection at any time 
during a one hour launch window. They also allow for variations in 
thrust, specific impulse and weight in any of the stages as well as 
perturbations due to winds during Stage 1 flight. /J j/TT/M 


INTRODUCTION 

An empirical method of implementing the path adaptive guid- 
ance scheme involves the development of steering polynomials. These 
polynomials must be able to steer the vehicle toward a predetermined 
end condition in an optimum manner. 

The purpose of this, study is to calculate a set of steering 
polynomials capable of acceptable injection into a .100 n.mi. orbit. 
These polynomials allow for launch at any time during a one-hour 
launch window. Perturbations due to thrust, specific impulse and 
weight variations in any of the stages and the effect of winds during 
first stage are included in the development of the polynomials. 

The assumptions involved in the study are as follows: 

1. Saturn V Vehicle 

The launch vehicle is a typical three- stage Saturn V. Only 
that portion of the third stage required to achieve the proper orbital 
conditions is utilized during boost. The remainder of the propellant 
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is retained for use during insertion into the lunar trajectory. 

Vehicle characteristics used in the study are given in the following 
table. 

Stage 1 Stage 2 Stage 3 

7.500.000. (s.l.) 1,000,000. (vac)200, 000. (vac) 

6 . 000 . 000 . 1 , 366 , 078 . 359 , 267 . 

(minimized) 


Thrust (lbs) 

Liftoff Weight (lbs) 
Propellant Weight (lbs) 
Specific Impulse (sec) 


4,224,210. 


919,011. 

• Classified — 


2. Trajectory Optimization 

All of the trajectories used to generate the polynomials 
have been optimized for maximum burnout weight. 

3. Spaced-fixed Waiting Orbit 

The target orbit is a 100 n.mi. circular orbit with an in- 
clination of 28.52 degrees. The space-fixed orbit is oriented so 
that its nodal line is coincident with the line of nodes of the lunar 
plane at the date of launch. 

4. Launch Window 

A launch window one hour in duration was selected for this 

study. 


ANALYSIS 

The analysis consists of two principal parts: 

1. Determination of optimum trajectories which meet the 
desired end conditions. 

2. Fitting pitch and yaw steering polynomials which steer 
the vehicle along these optimum trajectories. 

1. Determination of Optimum Trajectories 

A zero lift trajectory is flown during Stage 1 flight. 
However, during the early portion of the flight a pitch maneuver must 
be performed to turn the vehicle from the vertical. This pitch 
maneuver, as well as the launch azimuth, must be optimized. No thrust 
is applied in the yaw direction during Stage 1. Consequently, except 
for perturbations due to oblateness and earth* s rotation, the first 


8 


stage flight is two-dimensional. 

Upper stage trajectories are optimized by a three-dimen- 
sional calculus of variations analysis with an oblate earth model. 
The optimization criterion for all phases of the flight is maximum 
burnout weight. 

The desired cut-off condition corresponds to the level of 
energy and angular momentum associated with a 100 n.mi. circular 
orbit over a spherical earth. On an oblate earth, this will result 
in an orbit whose radius from the center of the earth will vary in 
a periodic manner. 

Discussion of the optimization of the first stage pitch 
maneuver involves several variables which require definition. 

These are as follows: 

a) Tilt angle (CXp) ^This is a measure of the rate at 
which the vehicle is turned away from vertical attitude. A pre- 
programmed turn is made early in first stage flight with fol- 

lowing the pattern indicated in the following sketch. 


o c 



The values for t^, t^, are ^ xec ** Conse- 

quently, the only parameter which affects the M steepness M of the first 
stage flight path and the ensuing first stage burnout conditions is 

<V 

b) Launch Azimuth (A )/vThe vehicle azimuth at launch is 
measured clockwise from north. For purposes of this study, it is as- 
sumed that the pitch. plane of the vehicle is always aligned with the 
launch azimuth. This eliminates the need for considering any roll 
during the vertical rise portion of the trajectories. 
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c) Co-nodal time At the instant of liftoff the launch 
azimuth and latitude define a hypothetical orbit plane. If an azi- 
muth of 90° is used, there is one time each day when the line of 
nodes of the instantaneously established plane will be coincident 
with the line of nodes of the target orbit. This time is defined 
as the co-nodal time. 

d) Launch time (At.WThis parameter is a measure of the 
location within the launch window of the launch site at the instant 
of liftoff. For this study the space fixed target orbit has an in- 
clination equal to the launch latitude. This assumption places the 
launch site in the plane of the target orbit only once each day. 

This time corresponds to the co-nodal time defined previously and 
establishes an unique point, i.e., A t L = 0. Early and late launches 

are respectively characterized by negative and positive values of 
A tL. In the following sketch, A tL is depicted as an angle. This 
is the angle through which the earth rotates between the time of 
vehicle liftoff and the co-nodal time. 
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e) Plumbline coordinate system ''-'This coordinate system has 
its origin at the center of the earth. The yp axis is parallel to the 
gravity gradient which passes through the launch site. The x p axis is 
parallel to the launch azimuth and the z p axis is perpendicular to the 
x p -y p plane forming a right hand system. 

f) Steering Angles (X iX »"X- ) Eulerian angles which 
locate the missile axis with respSct tS the r plumbline coordinate 
system. 


r 
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The launch window studied is one hour wide, i.e., -30 min^ 
A tL ^ + 30 min. For each selected value of Ati, the parameters 
A z and Ot p must be optimized. Studies indicated that the optimum 
value of otp is quite insensitive to variations in A z and A tj^. 
This was a fortuitous result since an actual first stage flight may 
use an open-loop X-p versus time guidance which is invariant with 
launch time and launch azimuth. A fifth order polynomial in time 
for X p was written to fit the nominal X p versus time table. This 
polynomial which was used to guide all of the first stage trajec- 
tories flown in this study is shown below. 


First Stage x. Polynomial 


? 

J> > .5 

Xp = a o + a it + a 2 t + 

a^t + a^t + a<-t 

a = +0.3835668 
O 

_4 

a 3 = -0.9412463 x 10 

a 1 = -0.2384704 

a^ = -0.6213386 x 10 -7 

a 2 = +0.1633698 x 10' 1 

O 

a c = +0.1104191 x 10" 
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For each value of At L studied, trajectories were run at 
several launch azimuths to determine the weight placed into orbit by 
a nominal vehicle. These trajectories followed the ")£_p polynomial 
during first stage and were optimized in the upper stages by the cal- 
culus of variations to obtain maximum burnout weight. 

A similar analysis in Reference 1 discusses a parametric 
study involving variations in launch time, launch azimuth, target 
orbit inclination and altitude on a spherical earth. 

Figure 1 illustrates the results of this study by plotting 
A z versus weight in orbit. For each Ati an optimum value of launch 
azimuth is clearly defined. The optimum launch azimuth associated 
with each value of At L is given in Figure 2. Note that the optimum 
azimuth for Ati = 0 is 89.15° and not 90° as one might expect. This 
result is due primarily to the fact that the powered trajectory re- 
quired a finite time (about 12 minutes) and that the analysis was 
performed using an oblate earth model. As shown in the following 
sketch, the gravity gradient is resolved into two components, one 
directed toward the center of the earth and one perpendicular to the 
geocentric radius. 
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A vehicle launched due east (A z = 90°) at A tL = 0 would 
tend to leave its plane because of the gravity component, gl . In 
order to maintain flight in the desired plane, thrust would oe re- 
quired in the yaw direction. In effect, this would be a three dimen- 
sional trajectory. When the launch direction is slightly north of 
east, gravity draws the vehicle into the target plane without any 
expenditure of propellant. Consequently, although a small plane 
change is performed, the total energy expended in reaching the de- 
sired end condition is a minimum. 

In Figure 3, the maximum weight which can be placed in 
orbit at each value of AtL is given. This curve indicates that the 
best launch time with respect to weight in orbit is & tjj = +17 min. 
This corresponds to a launch azimuth of 91.9° and appears contrary 
to the expected optimum of 90°. Explanation of this result involves 
a trade between the effect of the gravity component, g j , previously 
discussed and the contribution of the earth’s rotation zo the total 
energy imparted to the vehicle. As launch azimuth is varied from 
north to south, the effect of the gravity component becomes more 
pronounced. In fact, if this effect is the only one considered, a 
launch direction of due south is optimum. However, the tangential 
velocity imparted to the vehicle due to earth's rotation is a maximum 
when the launch direction is due east. When the effects of earth's 
rotation and gravity are considered simultaneously , the optimum 
launch azimuth becomes some angle south of east. Since the earth's 
rotational effect is so much larger than the gravity effect, the 
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optimum launch direction is slightly south of east. 

Empirical implementation of the path adaptive guidance 
scheme involves forcing the vehicle to follow an optimum path. The 
guidance system continually monitors the vehicle's state variables, 
i.e., position and velocity components, thrust acceleration, time, 
and The on-board computer carries prestored guidance polyno- 

mials which allow calculation of the optimum required thrust direc- 
tion as a function of the state variables. Prior to flight, a volume 
of optimum trajectories is generated using the calculus of variations. 
These trajectories, when flown, place the vehicle at the desired end 
condition in an optimum manner. Polynomials are written to fit the 
optimum "X- history of these trajectories. 


Two basic considerations were involved in the selection of 
optimum trajectories contained in the volume used to write the poly- 
nomials. These considerations are size and content. It is desirable 
to make the volume large enough to adequately cover the range of 
variables required and small enough to prevent the solution from 
becoming too cumbersome. It is also desirable to vary the vehicle 
and flight parameters to assure that any given vehicle, whose charac- 
teristics are within their allowable tolerances, will fly a trajec- 
tory enclosed in the volume studied. The range of dispersions used 
is as follows: 


Stage 1 

Thrust. + 3 % 


Stage 2 Stage 3 

+ 3 % + 3 % 


Isp 


+ 4 sec. 


+ 6.36 sec. 


+ 6.36 sec. 


Stage one trajectories were also perturbed by the presence 
of head, tail, and cross winds. The wind profile assumed is maximum 
design wind at an altitude of 12 Km. This wind is defined in Refer- 
ence 2 . 


The volume chosen included vehicles with off-nominal values 
of thrust and specific impulse in each stage, taken one at a time; 
i.e., when Stage 1 was assumed off-nominal, Stages 2 and 3 were flown 
as nominal. In addition, extreme variations were included by assum- 
ing that the vehicle parameters all varied simultaneously in the same 
direction. For instance, negative dispersions in both thrust and 
specific impulse in all three stages were considered. The trajec- 
tories flown with winds included vehicle variations in first stage 
with all nominal upper stages. For each vehicle variation studied, 
optimum trajectories are flown for -30 min ^ ^ tp, £ +30 min. in ten 
minute intervals. This volume resulted in a total of 339 optimum 
trajectories. 
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Table 1 illustrates the range of first stage burnout condi- 
tions generated by the volume of vehicle and flight dispersions. The 
increments of position and velocity vector components referenced to 
the completely nominal vehicle are shown. An indication of the large 
size of the volume may be obtained from the maximum and minimum vari- 
ations denoted by asterisks. 

2. Curve Fit 


Once the volume of optimum trajectories has been calculated, 
tuc next step involves generating guidance polynomials capable of 
steering the vehicle during boost to an acceptable orbit injection. 
Pitch and yaw steering polynomials must be written. In addition, a 
polynomial which calculates the required time remaining to burnout 
during the latter portion of the third stage flight is also written. 

Steering and cutoff polynomials for the Saturn C-l vehicle 
are presented in Reference 3. In this study, liftoff always occurs 
when the launch site is in the target plane . Three dimensional 
turning during upper stage flight is required primarily to eliminate 
out-of-plane perturbations caused by the earth's rotation and oblate- 
ness. 


The polynomials are of the following form: 
*x= A + A X + A Y + A Z + Ai X + A Y + . . . 

O 1 c. j *+ 

= f (x, y* z » x i y» z » f /m, t, a t L ) 

JL = f(x, y, z, X, y, z, F/M, t, A t L ) 

t c = f(x, y, z, x, y, z, F/M, t, A t L ) 


where x, y, z, 


x, y, z, 

F/M 

t 

At L 
p’ X y 
t 


components of position vector in plumbline coordinate 
system (meters) 

components of velocity vector in plumbline coordinate 
system (meters/sec) 

2 

thrust/mass (meters/sec ) 
time from liftoff (seconds)- 
launch window (minutes) 

Eulerian angles defining thrust vector direction in 
the pitch and yaw planes. (degrees) 

time remaining to engine cutoff (sec) 
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All combinations of the nine state variables to third order, plus a 
constant, result in 220 terms, 

Selected points from the 339 optimum trajectories in the 
volume are used to write the steering and cutoff polynomials, A 
cross-section of the 220 possible terms is chosen up to a maximum 
of 50. Using those terms, a least squares curve fit technique is 
utilized to generate the desired polynomials. 

RESULTS 

A number of different polynomials were calculated for pitch, 
yaw, and cutoff. Separate polynomials were calculated for second and 
third stage. A preliminary comparison between polynomials may be 
based on the root-mean-sum of the difference between as defined by 
the optimum trajectories and ^ as calculated by the polynomial. Past 
experience indicated that if RMS is less than approximately 0.3* the 
polynomial would be acceptable and would result in an acceptable tra- 
jectory. In the course of running actual trajectory simulations with 
guidance commands provided by the polynomials, it was found that RMS 
is not too valuable an indication of their validity. Comparison be- 
tween two polynomials often indicates that the one with the lower 
value of RMS results in greater deviations from the desired end con- 
dition. 


Typical polynomials are presented in Tables 2 through 6. 
The coefficients for each variable are given. The term Y found in 
the polynomials is defined as Yp - R where Yp = value for Y in the 
plumbline coordinate system and R = radius of earth. The number of 
terms in each polynomial and its corresponding RMS are tabulated 
below. 


No. of terms RMS 


Second Stage Pitch 

48 

.136 

Deg. 

Second Stage Yaw 

45 

.081 

Deg. 

Third Stage Pitch 

45 

.854 

Deg. 

Third Stage Yaw 

45 

.191 

Deg. 

Third Stage Cutoff 

58 

.118 

Sec. 


In order to adequately test the validity of the polynomials, 
a number of 3-dimensional digital flight simulations were run. These 
simulations actually used the steering and cutoff polynomials to 
guide the flight path. The check trajectories were chosen to encom- 
pass each of the variations used to generate the initial volume. 
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The desired cutoff conditions are as follows: 

Velocity (V) = 7794.58 m/sec. 

Radius (r) = 6555200. meters 

Flight Path Angle (0) = 90° 

Orbital Inclination (i) = 28,52° 

The results of this check, using the polynomials already 
presented, are given in Table 7- Errors in altitude, velocity, 
flight path angle and orbit inclination are presented for each 
vehicle dispersion. The RMS of each error is also given. 

An indication of the amount of plane change required for 
liftoff anywhere in the launch window is presented in Figure 4. The 
angle A <Q is defined as the angle between the target plane and the 
plane established if an orbit is instantaneously attained at any 
given launch latitude and azimuth. Data are given for launch azimuth 
from 86 degrees to 94 degrees and for the optimum launch azimuth as- 
sociated with each value of At^. 

Ground tracks of boost trajectories for A tg = 0 and 
-30 min. are illustrated in Figure 5. Eor comparison, the ground 
track of the target orbit is also shown. The tracks are plotted in 
Mercator projection on a non-rotating earth. 

CONCLUSIONS AND RECOMMENDATIONS 

1. When launch occurs at the co-nodal time, the optimum 
launch azimuth is 89.15° • Maximum weight is placed in orbit if 
launch is delayed until 17 minutes after the co-nodal time. 

2. Guidance polynomials have been written which success- 
fully steer the vehicle into a space fixed orbit. Liftoff may take 
place anywhere within a one hour launch window. 

3. Perturbations in the nominal boost trajectory due to 
vehicle variations in thrust, specific impulse and weight are easily 
handled by the polynomials. Perturbations due to winds during boost 
are also readily handled. 

4. Simulation of typical boost trajectories steered by 
the reported guidance polynomials results in RMS errors of approxi- 
mately 1 km. in altitude, 0.1° in flight path angle, 0.4 m/sec. in 
velocity, and 0.007° in orbit inclination. 

5. The selection of terms for use in the polynomials has 
involved a trial and error process. Errors considerably smaller than 
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those presented can result if the proper selection of terms is made. 
Development of an efficient scheme for determining the most signifi- 
cant terms is desirable. 

6. Curve fit techniques other than least squares should 
be investigated. 

7. Comparison between polynomials on an KMS basis has 
been found to be ineffective. Consequently, in this study comparison 
has been based only on a digital flight simulation using the poly- 
nomials to steer the vehicle. An attempt to find a better statisti- 
cal comparison should be made. If some polynomials can be eliminated 
by the use of a statistical comparison, the amount of digital compu- 
ter time required to determine the best combination of terms in the 
polynomial would be greatly reduced. 


TABLE 1 


H 

00 


FIRST STAGE CUT-OFF CONDITIONS 
At^ = -30 Min. A^ = 86.1 Deg. 


RUN 

NO. 

FIRST STAGE 

DISPERSION 



FIRST STAGE 

CUT-OFF 



F 

Isp 

Wind 

AX 

(m) 

AY 

(m) 

AZ 

(m) 

AX 

(m/s) 

AY 

(m/s) 

AZ 

(m/s) 

/S 

CsecJ 

A-l 

+ 3 % 

+ 1 . 3 % 

None 

- 1133 

+3971 

- 56** 

+20.93 

+70.60 

- .28 

- 2.1 

A-15 

+ 3 % 

- 1 . 3 % 

it 

-10398 

+1808 

- 167 

-85.65 

+81.65 

- .005** 

- 6.5 

A-43 

- 3 % 

+ 1 . 3 % 

m 

+11157 

-2311 

+ 177 

+85.29 

-86.13 

- .02 

+ 6.9 

A-57 

- 3 % 

- 1 . 3 % 

ii 

+ 1014** 

-4060 

+ 58 

-24.94 

-69.92** 

+ .08 

+ 2.3 

B-l 

+ 3 % 

+ 1 . 3 % 

Head 

- 2294 

+4365 

- 112 

+ 8.63** 

+74.39 

- .82 

- 2.1 

B-3 

+ 3 % 

- 1 . 3 % 

n 

-11517 

+2186 

- 221 

-98.28* 

+85.45 

- .55 

- 6.5 

B-34 

- 3 % 

+ 1 . 3 % 

Tail 

+ 12228 * 

-2704 

+ 227 

+95.74 

-89.22* 

+ .41 

+ 6.9 

B-36 

- 3 % 

- 1 . 3 % 

rr 

+ 2046 

-4442* 

+ 106 

-14.38 

-73.10 

+ .51 

+ 2.3 

B-57 

+ 3 % 

- 1 . 3 % 

Lf .Cross 

-10474 

+1807 

+1105 

-86.76 

+ 81.58 

+13.44 

- 6.5 

B-6l 

- 3 % 

+ 1 . 3 % 

ii ii 

+11072 

-2312 

+1588* 

+84.38 

-86.22 

+ 13.22 

+ 6.9 

B-82 

+ 3 % 

+ 1 . 3 % 

Rt .Cross 

- 1121 

+3946 

-1372 

+20.99 

+70.33 

- 13 . 62 * 

- 2.1 

B-84 

+ 3 % 

- 1 . 3 % 

ii ii 

-10386 

+1784** 

-1434 

- 85.82 

+81.38 

-13.41 

- 6.5 


Nominal Case 

x = 131216 

Y = 6435131 

Z = 22224 

* Maximum Variation 


x = 2563.69 

Y = 855.93 

z = 17.71 

** Minimum Variation 


*c = 153.7 

Xp = 56.6° 

X y = - 791 ° 



TABLE 2 


YAW STEERING POLYNOMIAL 
SECOND STAGE 


TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

A (Deg) 
o 

-9.3163288 

YT 

O 

-7.8257465x10 

XY F/M 

t5. 094291 3xlO“ 9 

X (M) 

+1.0195762x10~ 5 

ZT 

-2.7299943x10 -6 

YZ F/M 

+9.0348244x10“ 9 

Y (M) 

+3.1438830x10 -5 

Xt 

r 

+3.7748143x10“ 

z 2 f/m 

+5.l456394xl0“ 7 

Z (M) 

_A 

+8.2167335x10 

. YT 

f. 

-8.9242263x10 

y(f/m) 2 

O 

+1.1035070X10" 

X (M/S) 

-1.6089329xl0~ 3 

ZT 

-2.5124504x10“ 3 

z(f/m) 2 

0 

-4.1817956x10“ 

Y (M/S) 

-7.4887707xl0 -3 

F/M T 

-3.7529165x10"^ 

2 

Y x 

_l4 

- 4 . 0329277 x 10 

Z (M/S) 

+3.5143001x10 -3 

T 2 

+4. 5386163x10“ 5 

XZT 

+9.5624569x10“ 13 

F/M(M/S 2 

)+2.4306434x10 _1 

YAt L 

-7.3868989x10“ 7 

YZT 

+4.7792121x10“ 12 

T (Sec) 

+6.0046996x10 -3 

z^t L 

-7 . 1529175x10"^ 

z 2 t 

+2.1335529xlO“ 12 

t L ( Min )+3 • 283686lxlO _1 

tAt L 

-2. 6115389x1 0“ 3 

XZAt L 

+5 . 93358o4xlO -12 

YZ 

-1.0479217x10 -9 

ZAt L 

-6.7082010x10“ 5 

YZ At L 

+2.0924695X10 -11 

ZY 

+2.4622173xlO -7 

F/M At. 

-9.278l489xlO“ 4 

z 2 At L 

-1. 7755862x1 O -11 

*t 

+1.0077969x10"^ 

TAt L 

-6. 0591885x1 0 _if 

XY^t L 

+6 . 9803377x10 -11 

YZ 

-4 . 0146217x1 0~ 7 

XY F/M 

-1 .6504575x1 o“ 13 


+1. 3324114x1 0“ 10 

XT 

O 

-3.8593084x10 

YZ F/M 

-4 . 38l9295xlO -11 

ZY At T . 

-1.3949240x10“ 10 

Y = Y 

- R where R = 6373337*2 meters 
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TABLE 3 

PITCH STEERING POLYNOMIAL 


ro 

o 


SECOND STAGE 


TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

A 

-2.5197237xl0 +1 

YAt T 

+1.19l8824xlO " 6 

z 2 f/m 

- 1 . 05719 l 8 xl 0" 12 

o 

X 

+1.4428352x10”^ 


c 

+7.9003317x10 

YX F/M 

-8.4894779xl0 " 10 

Y 

+8.5079440x10 -5 

F/MAt T 

- 3 . 2709072 xl 0" 4 

Z(F/M ) 2 

-8.9403287x10 -7 

Z 

r 

+ 3 . 2921825 x 10 " 

TAt T 

+1.1509458xlO " 5 

p 

X^T 

-1.4671355x10" 13 

Y 

+1.5740623xl0 -2 

XY 2 

+ 1 . 6779022 x 10 

XYT 

+4.5312886x10 -19 

• 

z 

-4 

+3.1747471x10 

Y 3 

-I XT 

+4.6502831x10"' 

YYT 

+4.2224126x 10" x0 

F/M 

+7.7751782 

x 2 y 

-5.3540472x10" 15 

XYT 

+3.5908920x10" 9 

^ ^L 

+2.343392 / xlO "- 5 

YY 2 

+5.6684405x10 _11 

YZT 

+ 5 . 8405668xl0 -11 

Y 2 

+3.1233771xl0 " 10 

Y 5 

+ 5 . 0011765 xl 0" 9 

2 

XT 

+ 5 . 0728903 x 1 o " 10 

ZY 

+2 . 8067417x10"^ 

y 2 z 

+4.7123697xl0 " 14 

. 2 

XT 

+1.47962l7xl0 " 7 

YZ 

-3 . 3963735x10"^ 

z 2 z 

- 1 . 8458686x1 o " 1 2 

• 2 

YT 

- 1 . 8725866 xl 0" 7 

ZZ 

+3.4l87530xlO " 8 

• p 

YZ^ 

-1.5884291X10 " 11 

F/M T 2 

-4. 2467833x10 " 9 

z 2 

+4.3455238x10" 6 

*2 

ZZ 

+ 2 . 8702789 xl 0" 10 

t 3 

- 8 . 6560105 xl 0" 7 
- 1 

x at T 

- 3 . 3 l 88386 xlO " 9 

• • 2 

YZ 

+1.254238lxl0 " 9 

XYAt 

+2.0395705xl0 _i 


- 2 . 2090197 xl 0 -8 

XY F/M 

-7.l365936xl0 " 12 

2 L 
Y 2 At L 

-1 4 

+4.1764537x10 

S4t L 

-5.4797011x10" 8 

y 2 f/m 

-2.0949745xl0 " 12 

XXAt L 

+ 8 . 1906996 xl 0" 13 


Y = Y - R where R = 6373337.2 meters 092463 

P 


TABLE 4 

YAW STEERING POLYNOMIAL 
THIRD STAGE 


TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

A 

o 

-4.8682864 

YT 

O 

+2.6615327x10" 

XY F/M 

+2 . 5244878x10"^ 

X 

+1.8l70245xl0 -5 

ZT 

-3.727123^x10 

YZ F/M 

0 

+4.2969013x10" 

Y 

+1 .7919258X10" 3 

XT 

zr 

-3.7089316x10" 

z 2 f/m 

r 

+4.8997438x10" 

Z 

+1.352?021xl0" 3 

YT 

-1.5211l84xlO" 5 

y(f/m) 2 

+1.7803666X10" 7 

i 

+1 . 57 z +2798xlO -3 

ZT 

_4 

+7-5556900x10 

z(f/m) 2 

+4. 6359573x1 O" 6 

• 

Y 

+1.2294579x10 -2 

F/MT 

_4 

+2.4599282x10 

2 

YT 

-8.7370124x10" 15 

# 

z 

-3.7419080x10 -1 

T 2 

_4 

+1.1840582x10 

XZT 

+3. 3831325x1 O" 13 

F/M 

-9.0590695X10' 1 

YAt L 

-7.6782971xlO" 7 

YZT 

+4. 70317^5x1 O" 13 

T 

-4.3223290x10 -2 

ZAt L 

-9.6857025xlO" 7 

Z 2 T 

-4.9203751xl0" 13 

At L 

+2.3549313x10 -1 

YAt L 

-4 

-2.1159120x10 

XZAt L 

-1 4 

-7.7466600x10 

YZ 

-1.22l8951xlO" 9 

zat L 

+2.0292869xl0" 4 

YZAt L 

-1.8050l6lxl0" 12 

ZY 

0 

+4.5606823x10" 

F/M At 

+3. 91^6195x1 o" 2 

z 2 /u l 

-2. 9903261x1 O" 1 2 

• • 

XY 

-6 .0549310x1 O" 7 

TAt L 

-1.09723l6xl0" 3 

XYAt L 

+1. 1857600x1 o" ]0 

YZ 

+8.7691277xlO" 8 

XY F/M 

-5.85198 / i2x10" 13 

YYAt L 

+2. 185979^x1 O" 10 

XT 

-2.53l8894xlO -3 

YZ F/M 

-4. 2464884x1 O" 11 

ZY AtL 

+1. 8029512x1 O" 10 


Y = Y - R where R = 6373337 • 2 meters 
P 
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TABLE 5 

PITCH STEERING POLYNOMIAL 


ro 

ro 


THIRD Stage 


TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

1 

TERM 

COEFFICIENT 

A 

o 

-4.0546367x10 +5 

• • 

XY 

+1 . 9631275x10”^ 

XY F/M 

+1 .0527299x1 o ” 12 

X 

-1.9664022x10" 5 

Y 2 

-4 

+3.4629657x10 

YX F/M 

0 

-2.1355763x10” 

Y 

+7.394l658xlO - ^ 

Y F/M 

—4 

+1.8056825x10 

2 

XT 

+2 . 5336599x1 0 ” 12 

Z 

- 2 . 99628 l 3 xlO ' 5 

XT 

-2.9404846x10“ 6 

XYT 

+4.12l6l78xlO ” 12 

X 

+2.5675969X10" 1 

YT 

-3 . 3249322x1 0” 9 

2 

Y^T 

-7.1577965x10” 12 

• 

Y 

-4.9720504x10 _1 

XAt L 

+3.643889x10 -8 

*2 

XT 

0 

+7.3205544x10” 

• 

Z 

+2.6498565x10”^ 

YAt L 

0 

- 9 . 1311300 x 10 ” 

YYT 

+2.6467831xlO ” 10 

F/M 

-7.4008585 

ZAt L 

+ 1 . 0703024xl0 ~ 8 

XYT 

0 j 

- 7 . 3210298 x 10 " 

T 

+ 1 . 372 l 689 xlO +1 

YAt L 

+2. 295906 9xlO -9 

2 

XT 

-9.8l67372xlO ” 9 

At L 

+4. 9925804x1 0” 2 

zAt L 

- 1 . 9619938 x 1 o _if 

• 2 

XT 

-2.5073155xl0 ” 6 

XY 

+5.2647937xlO -9 

F/M At L 

+2.6763259x10“^ 

F/M T 2 

+4.4789852x10 -5 

Y 2 

0 

+1.0047664x10” 

TAt L 

.4 

-1.7694321x10 

t 3 

+I.606l42lxl0 ” 5 

XX 

+4.2551235xlO " 7 

x 2 5 

+ 6 . 700 l 801 xlO" lZt 

XYAt L 

-1.338583lxl0 ” 14 

YX 

-3.oi8o8i9xio” 8 

YY 2 

- 3 . 2219629 x 1 o ” 10 

Y >t L 

-l 4 

-5.8639033x10 

X 2 

-3.7H0768x10” 5 


+ 6 . 9654744x1 O -8 

xxAt L 

-1.9954698xlO -12 


Y = Y - R where R = 6373337 • 2 meters 092463 

P 


TABLE 6 

CUTOFF TIME POLYNOMIAL 
THIRD STAGE 


TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

TERM 

COEFFICIENT 

A 

o 

+2.4983B50x10 +4 

ZF/M 

-l .125604 2xl0 -3 

Z 3 

-3.43l46l6xlO” 9 

• 

X 

-1.0576426x10 +1 

(F/M) 2 

+2.4077407x10 +1 

X(F/M) 2 

- 3 .l 809302 xl 0” 3 

• 

Y 

-4.3829501 

XT 

+6.2334678x10" 3 

Y(F/M) 2 

+5.3836821X10" 4 

• 

Z 

+4 . 8057398xl0 _ ^ 

YT 

-4 

+2.7125073x10 

Z(F/M) 2 

+8.l422970xlO” 5 

F/M 

-4.9569286xio +2 

T 2 

-i . 3807829 x 1 o” 2 

*2 

XT 

-6. 4530036 x 1 o" 7 

T 

-1.3299996x10 +1 

XAt L 

-1 . 6312105x10”^ 

•p 

YT 

-9.7423446x10 -8 

At L 

-1.3622844x10 -3 

iAt L 

+2.1012l8lxl0” 6 

Z^T 

-3.69159^3x10” 9 

X 2 

+1 . 5076425x10”^ 

TAt L 

+1.4034l74xl0” 5 

• 2 

XT^ 

£ 

+2.1615007x10 

• • 

XY 

+1. 454201 3xl0” 3 

(At L ) 2 

_4 

+1.5633668x10 

• 2 

YT 

- 3 . 0226 l 62 xlO ” 7 

*2 

+2.5875335xl0” 4 

(X) 3 

O 

-7.0478323x10” 

• 2 

ZT 

+1.7156491x10” 10 

*2 

Z 

-5^8075911x10”^ 

x 2 y 

- 1 . 1930620 xl 0" 7 

•2 

Y 2 4t L 

+i.279i656xio” 10 

XF/M 

+6.2449482xio” 2 

XY 2 

O 

-3.6179801x10” 

(4t L )3 

+6.2766209x10 -7 

YF/M 

-1.7262958x10 -2 

Y 3 

-7.6967234xlO" 9 




Y = Y - R where R = 6373337.2 meters 
P 
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TABLE 7 


ro 

-f=- 


CUTOFF ACCURACY 

DESIRED CUTOFF CONDITIONS: Radius (r) = 6,555,200 Meters 

Flight Path Angle (•&) = 90.0 Degrees 
Velocity (V) = 7794-58 Meters/Second 
Inclination (i) = 28.52 Degrees 


RUN 

No. 

Xt L 

(Min) 

DISPERSION 

Ar 

(KM) 

A 

(Deg. ) 

A V 

(M/Sec. ) 

A i 

(Deg. ) 

1 

0 

Nominal 

+0.66 

-0.06 

+0.11 

-0.009 

2 

0 

-2% in Stage 1 Thrust 

+1.63 

-0.24 

+0.16 

-0.009 

5 

0 

-1% in Stage 1 Isp 

+1.05 

-0.08 

+0.34 

-0.005 

4 

0 

+10% in Stage 1 Inert Wt. 

+1.69 

-0.17 

+0.31 

-0.004 

5 

0 

-2% in Stage 2 Thrust 

-0.28 

+0.09 

+0.50 

-0.007 

6 

0 

-1% in Stage 2 Isp 

+1.48 

-0.14 

+0.20 

-0.009 

7 

0 

+10% in Stage 2 Inert Wt . 

-1.01 

+0.12 

+1.16 

-0.006 

8 

0 

-2% in Stage 3 Thrust 

+0.62 

-0.05 

+0.22 

-0.008 

9 

0 

-1% in Stage 3 Isp 

+0.67 

-0.06 

+0.11 

-0.009 

10 

0 

+10% in Stage 3 Inert Wt. 

+0.26 

-o.i8 

+0.39 

-0.007 

ll 

0 

Head Wind 

-0.55 

+0.06 

+0.52 

-0.002 

12 

0 

Left Cross 'Wind 

-0.82 

-0.07 

+0.13 

-0.008 

13 

-30 

A t T = -30 Minutes 

JLi 

+0.70 

-o.o6 

+0.11 

+0.002 

kMS 

0.99 

0.12 

0.43 

0.007 
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Summary 


A 


A set of guidance equations for power flight of the third stage of 
the Saturn V from parking orbit to earth- to-moon transit conic have been 
developed. This set consists of four polynomial expressions in terms of 
time and the state and vehicle performance variables. These equations 
provide the time of leaving the parking orbit, pitch and yaw steering 
angles during power flight and time remaining to third stage cutoff. 

The method utilized to obtain these guidance equations is based upon 
the calculus of variations. n u ~ # 
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INTRODUCTION 

The purpose and object of this report is to present the work in- 
volved and the results obtained to date concerning the development of a 
set of guidance equations for properly guiding a Saturn V third-stage 
vehicle to: (1) take off from a perturbed parking orbit, (2) boost 

through a 3-dimensional optimal twisted powered flight, and (3) inject 
into a properly selected earth-moon transit conic, so that the vehicle 
will free-fall into the vicinity of the moon, passing the moon at a pre- 
determined minimum altitude over a pre-determined point on the moon’s 
surface. These guidance equations would command the control system of 
the vehicle in terms of ’’when” and ”how” : (1) when to re-ignite the 

third-stage rocket, (2) how to steer the rocket thrust in pitch and 
yaw, and (3) when to finally cut off the rocket thrust, thereby, effec- 
tively carrying out the optimal 3-dimensional post-orbital boost to 
transit . 

Briefly, the report first describes and discusses: (1) the basic 

desired unperturbed optimal trajectory, (2) various unavoidable pertur- 
bations including analysis of perturbed parking orbits, and (3) the 
3-dimensional post-orbital boost which was studied following the per- 
turbed parking orbits and other perturbation studies. A large volume of 
perturbed trajectories was computed on the basis of all types of pertur- 
bation combinations. This spectrum of perturbed optimal trajectories 
was analyzed at every time interval and was reduced into a representa- 
tive statistical model. All controlling guidance parameters were ex- 
pressed in the form of multi-term polynomials of the vehicle’s state 
and performance variables. These guidance polynomials may be used to 
convert instantly sensed state and control elements of the vehicle 
into command signals to enable the vehicle to follow a newly selected 
optimal path from instant to instant. Finally, each of these polynomial 
guidance functions is analyzed: (1) in terms of its own accuracy as to 

how closely it represents all of the perturbed optimal trajectories by 
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the root-mean-square value of their value differences at all time inter- 
vals, and (2) in terms of its guidance accuracy as to how far off is 
the cutoff performance of the vehicle when it is steered by these guid- 
ance polynomials as compared with the cutoff performance if it is steered 
by the theoretical function of variational calculus optimization. 
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PRE-SELECTED UNPERTURBED NOMINAL OPTIMAL TRAJECTORY 

The desired and nominal trajectory chosen for this study represents 
a 72-hour space flight path to the moon on October 13, 1966. This date 

4 

was chosen because there was an angular difference between the minimum- 
inclination parking orbit from Cape Canaveral launch site (90° azimuth 
over Cape Canaveral) and the earth-moon plane, thereby illustrating the 
3-dimensional feature of the guidance functions to be developed in this 
paper. 

This nominal trajectory will pass the moon at a minimum radius ' 
(periselenum) of 1885 km (from the moon's center) over a predetermined 
point on the moon's surface, which is 5° latitude N. and 167° longitude 
E. from the zero point at the mean earth-moon- line. (Note: The said 

167° longitude E. is in the rear side of the moon as viewed from the 
earth; and the zero point at the mean earth-moon- line is in the Sinus 
of Medii which can be seen from the earth at all times.) By n-body 
celestial mechanics or by 6-body J.P.L. programming, a proper earth- 
moon transit may be approached by Cape Canaveral launching at lOh 36m 
6 sec (Greenwich time) on October 10, 1966. The vehicle will be in- 
serted eastward into a 90° azimuth 100-nautical mile parking orbit. It 
will then coast 1-2/3 orbit, and then inject into this earth-moon trans- 
sit for free flight to the moon, passing the pre-selected point in the 
vicinity of the moon. 

The point of interest involved in this project centered in the 
duration of the parking orbit waiting period and the post-orbital boost 
to earth-moon transit. The trajectory to be optimized will have both 
ends variable. It starts from the perturbed parking orbit in the 2/3- 
orbit region and ends on the time-varying earth-moon transit (see 
figure 2). The vehicle, under the adaptive guidance, will seek and 
follow an optimal path from instant-to-instant based on instaneous state 
and control variables sensed, until it reaches the transit conic and 
starts free flight. 

Figure 2 described this desired and preselected nominal and un- 
disturbed optimal trajectory. 
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PERTURBATIONS SURVEY AND PARKING ORBIT ANALYSIS 


It is usually impossible to fly along the above-mentioned pre- 
selected nominal optimal trajectory, because there are many unavoidable 
perturbations. To start with, there are perturbations due to a launch 
window of + 30 minutes for time to launch. Then, there are the initial 
errors of position and velocity vectors at the time of insertion into 
parking orbit. This insertion error, in terms of variations of velocity, 
altitude, and path angle, Av, Ar,A6> , will be modified and enlarged 
during the parking orbit coasting period due to earth oblateness and con- 
trolled venting effects of the vehicle's propellant tank. Therefore, 


parking-orbit coasting may be greater than at the time of insertion. 
Besides, there are perturbations due to the orbital boost window; i.e., 
one period advance or one period delay for the time of re-ignition. In 
other words, the vehicle may be ignited in the region of the 2/3 orbit, 
1-2/3 orbit, or 2-2/3 orbit. 

It is to be noted that the variation in yaw (out-of-plane motion) 
was assumed to be negligible for this first-cut investigation and was 
not included in parking orbit insertion errors. The out-of-plane yaw 
effect, however, was properly included in the studies of parking orbit 
perturbations due to earth oblateness and related factors. 

After the third stage re-ignition, there are further perturbations 
due to the imperfections of post-orbital-boost's propulsion system in 
terms of variations of thrust and specific impulse, +Af, + A I g ^. 
This, together with other perturbations as mentioned in previous para- 
graphs, will prevent the vehicle from following the exact preselected 
nominal optimal trajectory in reaching the final third stage cutoff 
point on the earth-moon transit. 


errors 


of Av, A r, at the time of re-ignition or at the end of 



The study of this complex picture of perturbations may best be 
approached by the study of the parking orbit perturbations. During 
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parking orbit coasting, there are two external variable forces acting 
on the vehicle. They are: (1) gravitational force due to the oblate 

earth, and (2) assumed axial thrust force due to the controlled venting 
of the vehicle f s propellant tank. As the pressure in the propellant 
tank rises, due to aerodynamic heating and solar heating, it must be 
released from time to time in order to prevent it from exceeding certain 
limits. Figure 3 illustrates the general venting mechanism of the 
third-stage vehicle; and figure 4 illustrates the force function due to 
the axial venting. The venting is assumed to start after 12 minutes 
from the time of insertion. It exerts a 300-pound thrust for 90 seconds, 
shuts off for 16 minutes, and vents again for 90 seconds. This venting 
cycle will be repeated until the time of re-ignition. 

The parking orbit perturbation due to initial insertion error was 
studied on the basis of this assumed venting-force function together 
with oblate earth effect. The plumbline coordinate system, the pitch 
plane diagram, as well as the orbital pitching rate of the vehicle 
adopted for the study, are illustrated in figures 4, 5, and 6. 

The initial insertion error in terms of Ar, Av, and Q was assum- 
ed to be + 1 km, + 2 m/sec, and + 0.02 degree. A total of ten trajector- 
ies for parking orbit coasting was studied to cover all cases of probable 
perturbation combinations with oblate earth and controlled venting ef- 
fects. One trajectory was for oblate effect only; one trajectory for 
combined oblate and venting effect; and the other eight trajectories 


and controlled venting vehicle. The details of the trajectory composi- 
tion for the parking orbit perturbation are listed in figure 7. 

Results of all ten cases were studied, analyzed, and compared. 
Evaluated results indicated that the initial insertion velocity error, 
Av, was most influential and represented the controlling input. Figure 
8 is a composite diagram designed to illustrate this point. Outer and 
inner ellipses (see figure 8) represent the energy range of the perturbed 


for various 



effects on the basis of oblate earth 
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orbit. They are based on initial ISv variations and they envelop 
other ellipses which are based on A r variations. With this controlling 
initial A v at insertion as input, the resulting perturbations of t v, 

A r , and A 0 in the re-ignition region were studied for a vented orbit- 
ing trajectory around an oblate earth. They are plotted in figure 9 
showing Ar, Av, and A6> deviations at the time of re-ignition of 
the third-stage vehicle due to the initial insertion error of Av 
(for a vented orbiting trajectory and an oblate earth). As shown in 
figure 9, an insertion Av error of + 2 m/sec will cause a total varia- 
tion of approximately + 10 km in altitude, + 8 m/sec in velocity, and 
+ 0.05 degree in path angle in the area of the third orbit re-ignition. 

In addition to this parking orbit perturbation (which may be repre- 
sented by velocity variations, + Av, at insertion), there are three 
other types of perturbations: launch window or launch time variation; 

orbital residence or waiting period variation-; and thrust I variation 

bp 

during third stage post-orbital boost. These are illustrated in figure 
10 and tabulated in figure 11 (table). This 3 set means that a total 
of 27 trajectories will give a fair coverage of all possible perturba- 
tion combinations for studies prior and up to third-stage re-ignition. 

In order to cover fairly the various perturbation combinations for 
studies of post-orbital boost, 27 x 5 = 135 trajectories are used. 

These perturbed trajectories may be seen pictorially in figure 10. 
Illustrated in figure 10 are launch windows, parking orbit windows, 
and Ar, Av, A0 at re^ignition due to insertion A v error. Also 
illustrated in figure 10 are the 3-dimensional post-orbital boost 
perturbations and the spread of the end earth-moon transit which is 
due to the combined effect of these perturbations and the shifting of 
the transit following the moon f s motion. 
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3-DIMENSIONAL TWISTED POST- ORBITAL- BOOST AND RE- IGNITION TIMING 

As previously mentioned, the post-orbital-boost trajectory may be 
optimized by variational calculus with both ends variable. It starts 
from the perturbed parking orbit and ends on the earth-to-moon transit 
trajectory at third stage final cut off. Such a space or osculating 
conic for nominal flight may be based on the nominal earth-to-moon tran- 
sit trajectory computed by the 6-body JPL program. To cover the 4-hour 
window of third-stage re-ignition (2/3 orbit, 1-2/3 orbit, 2-2/3 orbit 
with variations of early, nominal, and late launching for each orbit), 
a total of 9 earth-to-moon transit trajectories corresponding to differ- 
ent time of third-stage final cutoff and different time of arrival at 
Moon were computed and plotted in figure 10-A. The time-varying earth- 
to-moon transit may be further transformed into time-varying space or 
osculating conic passing the final cutoff point as illustrated in 
figure 10-B. The orbital elements of this osculating conic is tabulated 
in figure 10-C. 

As this table is further expanded to cover more osculating conics 
at smaller time intervals, these orbital elements may be plotted and 
curve-fitted as time varying functions. These functions, therefore, 
represent one of the two end conditions for the above-mentioned varia- 
tional calculus optimization. The other end condition, of course, may 
be represented by the perturbed parking orbit in the neighborhood of 
third- stage re- ignition. 

In view of the fact that the earth-moon transit plane forms only 
a very small angle with the due-east parking orbit plane from Cape 
Canaveral, the analysis for third-stage re-ignition may be approached 
by a 2-dimensional approximation. (See figure 12.) Through variational 
calculus optimization, both the time duration of the third-stage re- 
burning, t^, and the swept angle by the post-orbital boost > were 

computed on the basis of various cutoff points on the transit or various 
range angles between "cutoff points on the transit and perigee of the 
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transit, M (p . Graphs of t^ vs 
figure 13, indicate a first-cut 
ignition point. 


• CZ) and /j> . n vs. CD , as shown in 
/ p / inci / p 

information concerning optimal re- 


By this approximated location of re-ignition point, a rescheduling 
by the venting may be made. In order to provide a clear and smooth 
operation for the astronaut within 500 seconds prior to re-ignition, 
it was programmed to reschedule or to advance the last venting in 
order to have a 500-second no-venting coasting prior to the third- 
stage re-ignition. It is important, therefore, to record the space 
angular orientation of the vehicle at this particular time: time of 

"start of no-venting program" or t sonV p> the details of which are shown 
in figure 14. 


Simultaneously, the parking orbit perturbation was re-studied on 

the basis of this rescheduling of the last venting. Av, A r , and Ad 

variations at the time of "start of no-venting program" (t ) due 

A sonvp 

to initial insertion error of /Av were re-calculated. 


With the values of real time of t together with the corres- 

ponding Av, A r, and A U , several 3-dimensional calculus of varia- 
tions UOV) optimal paths were studied with different no-venting 
costing, At, along the last 500 seconds, plus or minus a few seconds, 
of the parking orbit. As also illustrated in figure 14, the optimal 

t ... was located after three such trials for each of the 27 
re-ignition 

trajectories to cover all probable perturbation combinations prior to 
re-ignition, as mentioned in previous sections and tabulated in figure 
11 . 


In detail, the 3-dimensional MSFC V-30 deck and calculus of varia- 
tions optimal program were used with state variables at real time t gQnv p 
and corresponding Av, Ar, and A $ as inputs. It started with an 
assumed parking orbit no-vent coasting of about 500 seconds, then pro- 
ceeded through re-ignition, and 3-dimensional optimal boost to final 
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cut off on transit. The procedures were repeated until the optimal 

t . . . was found for each of the above mentioned 27 trajectories 

re-ignition J 

(minimum time duration for post-orbital boost). 
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DEVELOPMENT OF GUIDANCE EQUATIONS: t . , X , X , t . ... 

x re-ig p y cut off* 

Once all 27 optimal trajectories are established to cover all prob- 
able combinations of perturbations prior to third-stage re-ignition, each 
of these trajectories may be examined prior to re-ignition at each of its 

time intervals. As the programed no-vent coasting from t lasted 

500 sonv P 

about 500 seconds, at 5-second intervals, a total of — x 27 or 2700 
points existed. At each of these 2700 points, we may note and tabulate 
corresponding time remaining to re-ignition together with vehicle’s 
state variables x, y, z and x, y, z. This is the guidance function of 
’’time remaining to re-ignition” expressed in tabular form. In order to 
handle the case by on-board computers, this large table is further re- 
duced by functional approximation through . least square procedure into 
32-term polynomials of state varialbes. In other words, we have de- 
veloped a statistical and representative model of the 27-element volume 
trajectory. In theory, this 32-term polynomial should fairly represent 
all of these 2700 points and should be capable of computing the "time 

remaining to third-stage re-ignition," t . , based on instantaneous 

r e — lg 

sensing of the vehicle's state variables. Figure 15 shows this 32-term 
equation. The units for x, y, z; x, y, z; and t are m, m/sec, and sec, 
respectively. 

An overall re-examination at each of these 2700 points was made 
in terms of the difference between the actual value and the value com- 
puted by the polynomial. This examination revealed that the root-mean- 
square value of all of these differences is equal to 1.3758 seconds. In 
a way, this RMS value indicates how well the polynomial represents these 
27-element volume trajectories and the corresponding 2700 points as 
their representative statistical model. 

Once this polynomial guidance equation for "time remaining to re- 
ignition" is obtained, each of the 27 trajectories may be re-run from 

the time of "t " through the polynomial-calculated time of re- 

sonvp 


44 


ignition, then, through the 5-perturbation combination for post-orbital 
boost (see figure 11). A total of 27 x 5 or 135 trajectories was com- 
puted. 

* 

Using the same procedure as described in previous sections and based 

on these 135 trajectories for post-orbital boost period of about 350 

350 

seconds at 5-second intervals, a total of approximately 135 x — — = 9000 
points may be tabulated. At each point, there are corresponding X^, X^, 
and M time remaining to final cut off," t and corresponding state 

and performance variables of the vehicle: x, y, z, x, y, z; and F/m, 

m/m, and t. By least-square curve fitting, a 43-term polynomial for X^, 
a 42-term polynomial for X^, and a 42-term polynomial for t^^ were 
developed and are tabulated in figures 16, 17, and 18. RMS values re- 
presenting the polynomials' accuracy are 0.455°, 0.105°, and 0.600 sec- 
onds, respectively. The units used for these guidance equations are 
2 — 

m, m/sec, m/sec , 1/sec, and sec for x, x, F/m, m/m, and t, respectively. 


* Note: Xp, X^ represents pitch steering angle with respect 

plumbline vertical and yaw steering angle measured from X, 
plane as shown in figure 6. 


to 

Y 
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ERROR ANALYSIS AND EVALUATION 

Finally, the accuracy of these guidance polynomials in terms of 
the difference between the cutoff position and velocity vectors of the 
vehicle when steered by these polynomials and those obtained if steered 
by theoretical functions of variational calculus optimization have been 
examined. 


In order to do so, 15 simulated flights were made, 5 from each 

orbit with 2-sigma value deviations of /\F , /\l , Aw, /\t, and zero deria- 

Sp 

tions (nominal). These 15 simulated trajectories, as tabulated in 
figure 19, were run by both adaptive guidance polynomials and by 
theoretical calculus of variations (COV) optimization. Each of the 15 
cases was examined in terns of cut-off r, v, @, and N (altitude, 
velocity, path angle, and angular difference between normals of trajec- 
tory planes or transit planes). These 15 values of differences were 
again expressed in RMS errors. Figure 20 illustrates the details of the 
comparison. Respective RMS errors at cutoff are about 2 km for altitude, 
5 m/sec for velocity, 0.06 degree for path angle, and 0.008 degree for 
transit plane deviation. 
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CLOSING REMARKS 


1. This study represents only the first cut investigation to 
illustrate the feasibility of the adaptive guidance technique in guid- 
ing the Saturn V third-stage vehicle through post-orbital boost to earth- 
moon transit. No attempt was made to carry out overall optimization of 
the problem. 

2. Attempts were made to illustrate the 3-dimensional feature of 
the guidance technique. There was angular difference between the park- 
ing orbit plane and the transit plane for each of the 135 trajectories 
studied. The plane difference ranged from 0.9 degrees up to 1.4 degrees 
with 1.2 degrees for nominal case. It is to be noted that none of the 
five reference planes; MEP (earth-moon plane), equational, ecliptic, 
transit, and parking orbital; are co-planar. 

3. It is to be repeated that this paper was based on: (1) assumed 

initial errors at parking orbit insertion, (2) assumed relationships for 
space conics at the third stage final cut-off, (3) approximated re- 
ignition timing. It is planned that the study will be repeated with re- 
fined and actual inputs in three areas: (1) actual errors instead of 
assumed errors at the time of parking orbit insertion , (2) ref ined 
time-varying functions of the corresponding space conic elements in 

the region of third-stage final cutoff , and (3) refined and iterated 
re- ignition timing . Area (1) is presently pursued by launching phase 
study group. Area (2) is to be pursued in a manner as illustrated by 
figures 10-A, 10-B, and 10-C. Area (3) will be pursued in a manner 
as described in the second half of the section n 3-Dimensiona 1 Twisted 
Post-Orbital Boost and Re-ignition Timing.” The optimal re-ignition 
timing will be determined by variational calculus with both ends 
variable. As mentioned previously, one end will be represented by 
last section of the refined parking orbit with adjusted venting sche- 
dule together with 500-second no-venting costing. The other end will 
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be represented by the refined time-varying-functions of the earth-to- 
moon transit conics. 

4. In error analysis, this accuracy was based on cutoff errors 
at final cut off on the earth-moon transit: 2 km in altitude, 5 m/sec 

in velocity, 0.06 degree in path angle, and 0.008 degree in trajectory 
plane deviation. It is expected that these errors will be further re- 
duced in the projected refined study. Furthermore, mid-course correc- 
tion will further correct whatever the final error may be in order to 
arrive at the moon in the desired manner. 
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TRAJECTORY COMPOSITION FOR PARKING ORBIT PERTURBATION ANALYSIS DUE TO INSERTION & OTHER EFFECTS 

(Listing of Initial Conditions) 




Case 

R(KM) 

V(M/Sec.) 

0 (Deg) 

Variation 

Venting 

1 

6558.4716 

7794.72 

90 


No 

2 

6558.4716 

7794.72 

90 

— 

Yes 

3 

6557.4716 

7794.72 

90 

-AR 

Yes 

4 

6559.4716 

7794.72 

90 

+AR 

Yes 

5 

6558.4716 

7792.72 

90 

-AV 

Yes 

6 

6558.4716 

7796.72 

90 

+AV 

Yes 

7 

6558.4716 

7794.72 

89.98 

-A0 

Yes 

8 

6558.4716 

7794.72 

90.02 

+A0 

Yes 





, AR . AV A0 


9 

6558.9716 

7795.72 

89.99 


Yes 





AR _ AV + A0 


10 

6557.9716 

7793.72 

90.01 

2 ~ 2 2 

Yes 



AR = 

1 Kilometer 





II 

> 

<3 

2 Meters/Sec. 





A0 = 

0.02 Degrees 
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ORBITAL ELEMENTS OF OSCULATING CONIC AT THIRD STAGE FINAL CUT OFF 


o 


Third 

Coast stage 

in final 


Parking cut off 

Orbit (Seconds) rp (km« ) 


E 

2/3 

orbit 

40785.997 

6603.1811 

S 

2/3 

orbit 

42580.437 

6601.5446 

L 

2/3 

orbit 

44205.445 

6601.5940 

E 

1-2/3 

orbit 

46198.599 

6602.9947 

S 

1-2/3 

orbit 

47990.947 

6601.3528 

L 

1-2/3 

orbit 

49617.133 

6601.4055 

E 

2-2/3 

orbit 

51611.055 

6602.8070 

S 

2-2/3 

orbit 

53400.350 

6601.1611 

L 

2-2/3 

orbit 

55028.666 

6601.2160 


Note 1: 
Note 2: 


a (km # ) e 

i 

(Degrees) 

216,352.78 .96947956 

28.587781 

216,364.77 .96948881 

28.416395 

216,361.80 .96948817 

28.600123 

216,370.05 .96948286 

28.587475 

216,283.78 .96947828 

28.416078 

216,378.48 .96949139 

28.599867 

216,385.56 .96948591 

28.587154 

216,408.32 .96949673 

28.415755 

216,403.23 .96949576 

28.599537 


E represents early launching; S, 
Times are measured from midnight 


Time of 

ft oj Perigee Passage 

(Degrees) (Degrees) hr« min, sec. 


8.5506383 

9.8240908 

11 

17 

15.394 

8.5362614 

10.151603 

11 

47 

9.676 

8.5814669 

10.398030 

12 

14 

14.743 

8.5958230 

10.737073 

12 

47 

27.997 

8.5767595 

11.083926 

13 

17 

20.188 

8.6271262 

11.310609 

13 

44 

26.431 

8.6411123 

11.650305 

14 

17 

40.454 

8.6152130 

11.986828 

14 

47 

29.592 

8.6728923 

12.223285 

15 

14 

37.965 


standard launching; L, late launching 
(Greenwich Standard) 
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LISTING OF SIMULATED TRAJECTORIES UNDER 

VARIOUS CONDITIONS OF PERTURBATION COMBINATIONS 
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*&£ 0.5 deg. [ 

( B ) NOMINAL ) 

PARKING ORBIT f 1 
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*T 1 0 
lavhchhA -30 min\ 
WINDOW v J 


3 


RESULTANT CONTROLLING 
PERTURBATION IN TERMS 
OF x,yz &*$,z AT 3rd 
STAGE RE-IGNITION AREA DUE 
TO AV AT INSERT I ON 



PERTURBATION DUE TO PARKING 
ORBIT RESIDENCE DURATION - 


(C) 

THRUST PERTURBATION 
DURING POST ORBITAL BOOST 


j orbit \ 

*T 4 ljo>-b/t>3 
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COEFFICIENTS OF 32 TERM POLYNOMIAL 
FOR TIME REMAINING TO 3rd STAGE RE-IGNITION 


PARAMETER COEFFICIENT 


CONSTANT 

- .13704652*10* 

X 

-. 339/3094*10 '• 

y 

-.9504-5607*10’' 

z 

*.1/532524 

X 

+ .74598325*/0 +t 

9 

- .221 10006 *10+ * 

z 

- .34006194*10* 

T 

*.70403094 

X* 

* .535867/8 *!0~' 

y % 

- . 92944665*10 '* 

z* 

+ . 02790/05*10 "• 

yx 

+ .15053849*10’+ 

zx 

* .69324013*10* 

X* 

- .7 9131437*10* 

yy 

* . 12401051*10 ~ 4 

Z Y 

- . 14744044*10 ’+ 


RMS Error 


PARAMETER COEFFICIENT 


y z 

* .44265428*IO l 

z* 

* . 24678161*10 + ' 

xT 

* .11220720*10 * 

yT 

- .11698904*10 3 

zT 

* .1214604) p/O' 1 

T z 

- .12930040*10 s 

X* 

- .25057946*10 “ 

y* 

- . 47954404*10 /7 

x*T 

- .16/9581/ * IO' 13 

z*T 

- .13650941*10-'° 

x*T 

75418146 */0’ $ 

y*r 

r. 8 570/ 534 */0’ 9 

z*T 

- . 37410572*10 

XT* 

- .15337463*10* 

9 T \ 

+ . 5/702550*10 

zT* 

- .91033639 */0’ 7 


Units - meters, seconds 

sec. 

Y transformation i 


Y=Yp - 8373337.Z 


VJl 


Figure 15 




COEFFICIENTS FOR 43 TERM POLYNOMIAL 
FOR X P STEERING ANGLE FOR POST ORBITAL BOOST 

(BETWEEN 3rJ STAGE RE-IGNITION AND FINAL CUT-OFF) 


PARAMETER 

COEFFICIENT 

PARAMETER 

COEFFICIENT 

CONSTANT 

+ .17228688*10° 

xy T 

- .20522858*10-° 

X 

+ .+2059390*10-* 

y±r 

+ .10565799* I0’ 9 

z 

+ .87211159*10-* 

xT z 

+ .59936783*10'" 

y 

+ .3299792! *10* 

xT z 

+ . 11032575 *J0~" 

z 

+ .26336817*10* 

yr z , 

+ .++5IZ29l*IO' ,e 

tin 

-.l65+8+l+*IO* 

(r/»)T* 

+ .8336+886*10-° 

y* 

-.80986205*10*° 

T 3 

- .+12396+9*10" 

zy 

-.9579 0635 *\0~ 7 

xzz 

-.12923623*10 ” 

yz 

-.5758+052* /O' 9 

z*z 

-.75/62203*10'* 

±* 

-.10765835*10'° 

xyz 

-.960+2316 *I0~ 9 

Y* 

-.+5905236*10'" 

yz (f/m) 

-. 8583+1 Z8*I0' 7 

yyt 

- .1002036 7*10'" 

y (F/y») z 

- .25781+38*10-° 

y* 

-.12+107+9 */0' 9 

zyr 

+ .92378919* 10** 

yz* 

-.1095333+* /O'” 

zT 

90507301 */0'* 

xy (F/») 

-.22+59561* 10** 

xyx 

+ .1+758270*10*° 

y*(F/») 

-.33685 72+ *10*° 

x*y 

+ .60773933*10*° 

xy (F/m) 

+ .21107992*10-° 

*y T 

+ .66852788*10*° 

yz(F/m) 

+ .87729871*10-° 

yT z 

- .322350+9*10** 

Z (F/m) z 

+ .21137200 */0‘ 3 

yy(F/yn) 

-.32+765/5*10-° 

X*T 

+ .16025891*10*° 

xyy 

-.92053916*10” 

xyT 

+ .++030386*10” 

yV 

+ ./ 6/06/65 *ZO 

yyr 

+ .9369 7273*10” 

RMS Error 

=.+ 5+48 / deg. 


Figure 16 


COEFFICIENTS FOR 42 TERM POLYNOMIAL 
FOR X y STEERING ANGLE FOR POST ORBITAL BOOST 

j or.i/^iifTMAl iaia triklAI 


PARAMETER 

COEFFICIENT 

PARAMETER 

COEFFICIENT 

CONSTANT 

-.55+5/576*10' 

x z (r/») 

* .3/ 3734/ 6* /O' 7 

X 

- .Z+569037*/0~ 4 

y z (F/m) 

-.2/37+055* /O' 7 

y 

8/7/7797 */0~* 

x z (F/m) 

- .13906632*10-* 

z 

~.l+229007*/0' 3 

x (F/tp) x 

-.19652365 * /O m * 

X 

* .1/7/1731 *IO 

y <F/*P 

* .12565766*10’* 

• 

y 

z 

- .1+326+70*10'' 

z (F/m) A 

- .32572526*10'* 

- .59530/5/ x /O' 3 

x (F/m) x 

- .360 5 9335* /O’ 4 * 

F/» 

-.73/70/73x10' 


*.10337563* /O' 4 

T 

-.ll775353*/0’' 

i (F/m) z 

* .9 150 9227* /O' 4 

XT 

*.7/+97336*/0' 9 

x z (F/m) 

-. +1717 6I3*I0- ,Z 

yT 

-.5206599S*i0~ s 

x z T 

+ .67472623*10 * 

xT 

*./6+6+223*/0' 8 

xyT 

-./ 3662406* /O-'* 
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COEFFICIENTS FOR 42 TERM POLYNOMIAL 
FOR TIME REMAINING TO 3rd STAGE FINAL CUT-OFF 
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TRAJECTORY COMPOSITION FOR ERROR 
ANALYSIS AT 3yd STAGE FINAL CUT-OFF 
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ERROR ANALYSIS 

AT 

3rd STAGE FINAL CUT-OFF 
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SUMMARY 



The minimum and maximum range capability of the S-IVB Stage 
and Service Module used as vehicles for powered abort from an expected 
volume of Saturn V boost trajectories has been determined. Abort is 
assumed to occur between the time of second stage ignition and a time 
near final lunar injection. Abort trajectories terminate at a re- 
entry altitude of 120 km with flight path angles of 94° and 99.5°. 


INTRODUCTION 


A U~r-i 4 


During flight of the Saturn V vehicle, from lift-off to 
lunar injection, malfunctions may occur which would necessitate an 
abort. Abort may be defined as the reaction to malfunction which 
requires the immediate abandonment of primary and secondary missions. 


are : 


Three flight modes may be considered for the Saturn V. They 


1. Normal flight - no malfunctions 

2. Abnormal flight - minor malfunction 

a. Primary mission may still be accomplished. 

b. Secondary missions may be accomplished if the 
primary mission cannot be completed. 

3. Aborted flight - serious malfunction 

a. Unpowered abort. 

b. Powered abort. 
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Using the Apollo vehicle, unpowered abort can be accom- 
plished with acceptable re-entry velocity and flight path angle for 
most sub-orbital abort initiation times. However, the low maneuvera- 
bility of the Apollo does not allow a great degree of flexibility in 
landing site selection. 

This study is concerned with powered abort using the S-IVB 
stage or Service Module. Initially, the study was directed toward 
obtaining a unique landing site accessible from all abort points. It 
soon became apparent that such a single site did not exist and the 
problem expanded to determining the range capabilities of the selected 
powered abort stages so that the minimum number of abort sites could 
be chosen. For a select group of landing sites, the final phase of 
the study will be to determine Path Adaptive Guidance steering poly- 
nomials for use during powered abort trajectories. 

It is the purpose of this abort study to establish the two 
dimensional range capability of the S-IVB stage and the Service 
Module (SM) when used as abort vehicles between the time of S-II stage 
ignition and injection into lunar orbit. The ground rules imposed 
during this study were: 

1) No coast would be considered from ignition of the abort 
vehicle until the re-entry point was reached. 

2 ) No cross-range capability would be examined. The study 

would be two-dimensional. 

/ 

3 ) Abort trajectories would originate from an envelope of 
likely boost trajectories as defined later in this report. 

4) Abort trajectories would terminate with acceptable re- 
entry conditions at an altitude of 120 km. 

The data contained in this report represents an extension of the 
study reported in Reference 1. 


DEFINITIONS 


h - altitude 

v - velocity 

t* - flight path angle measured clockwise from local 
vertical 

0 - range angle measured from launch to the vehicle re 

entry point along a spherical earth surface. 

OC - angle of attack 

^ - time derivative of oc 

time derivative of cJc 
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ANALYSIS 


The Powered Abort Study consists of two parts, each charact- 
erized by the vehicle to be used during powered abort flight. Part 
one used the Service Module as the abort vehicle and part two used the 
S-IVB stage as the abort vehicle. In order to begin the abort study, 
a boost trajectory envelope was defined. Eighteen boost trajectories 
were generated by varying thrust, Isp and stage inert weights of all 
three stages of the Saturn V vehicle. The variations used are tabu- 
lated below. 


Thrust (lbs) 

Isp (sec) 

Inert Wt. (lbs) 


Stage Variations 
S-1C 

+ 450,000 
+ 4 

+ 15,000 


for Boost Volume 
S-II 

+ 100,000 

+ 4 

+ 5,000 


S-IVB 
+ 2 0,000 
+ 4 

+ 2,000 


The descriptive characteristics of the base line Saturn V 
vehicle used in this study are given in Table I. 


The two boost trajectories selected, from which abort was 
initiated, were the two which most closely represented the boundaries 
on an altitude-velocity graph of the eighteen boost trajectories. 
These two boost trajectories form the boost trajectory envelope shown 
in Figures 1 and 2, and are termed "boundary” trajectories. 


For the Service Module, seven abort times were initially 
selected. Three of the abort times selected were during sub-orbital 
flight and four of the abort times were during super-orbital flight. 
Four sub-orbital abort times between second stage ignition and third 
stage ignition were selected for the S-IVB stage. 


The abort times selected were as follows: 


Service Module 

S-IVB 

167 

sec. 

167 

sec. 



275 

sec • 

395 

sec. 

395 

sec. 



503 

sec. 

650 

sec . 



790 

sec. 



834 

sec. 



874 

sec. 



946 

sec. 




Criteria 

Near S-II ignition 


Near S-II burnout 


Near S-IVB re-ignition 


Near S-IVB final cut-off 
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As the study progressed, it was discovered that, under the 
established ground rules, abort with the SM at 9^6 seconds could not 
be accomplished. A time of 906 sec. appeared to be near the limit of 
SM superorbital abort time and, therefore, was used in the SM super- 
orbital portion of this study. Later another time of 915 sec. was 
also studied in attempting to precisely locate the time limit. From 
each of the boundary trajectories, at each of the chosen abort times, 
it is desirable to fly the vehicle to a set of state variables as it 
approaches the earth's atmosphere. The set of variables is defined 
as a "re-entry window". These re-entry conditions define the terminal 
point of powered abort flight. For this study, the re-entry window is 
defined by an altitude of 120 km. and a flight path angle variation at 
this altitude from approximately 90° to 105°, depending upon the re- 
entry velocity. Abort trajectories were run on a 2-degree of freedom 
Calculus of Variations digital program which includes a mathematical 
scheme capable of isolating on pre-set conditions at burnout, (ref- 
erence 2) 


In this study the end conditions were initially selected to 
be a re-entry altitude of 120 km and a flight path angle of 9^ degrees. 
However, for the sub— orbital SM cases, isolations on 9^ could not 
always be attained. This inability of the SM to attain a 9^° re- 
entry from all abort points is due to the low thrust to weight ratio 
which is unable to modify the steep, high gpeed near ballistic path 
of the vehicle. Inability to isolate on 7# = 9 L \° for sub-orbital SM 
cases suggested that isolations from the entire boost trajectory at 
some higher ZA might be possible. A value of 99-5 degrees was selec- 
ted so that re-entry conditions for each powered abort flight were 
h = 120 km and lA = 9b or 99-5 degrees. 

The technique of obtaining the maximum and minimum range 
angles for each abort time and from each trajectory consists of sever- 
al steps: 


1. A series of explicit Calculus of Variation (C0V) tra- 
jectories were generated by specifying oCo an< ^ ®^o an< ^ using the 
transversality condition to specify oC o* This entire set of runs 
was made with velocity extremalized and with range angle an open 
variable, i.e., no explicit control or optimization of 0 was per- 
formed. 


2. Each of the above series of trajectories was examined 
throughout its burn time for points that approached the desired re- 
entry window c'onditions. This selected burn time and the correspond- 
ing ©Co and ^C 0 f° r that trajectory were then used as initial condi- 
tions for an isolation run on the exact desired re-entry window con- 
ditions. An "isolation" case is a series of explicit trajectory runs 
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in which the digital program automatically varies oC 0 and oc o until 
altitude at burnout equals 120 km and the flight path angle equals 
either 94° or 99*5° * whichever is specified. After an isolation was 
obtained, the burn time was changed to either a larger or smaller 
value and a new isolation of the re-entry conditions was attempted. 
This time incrementation procedure was continued to establish the 
minimum and maximum burn times which would yield isolation on the 
desired end conditions, for the S-1VB stage, the maximum burn time 
to produce desired re-entry was taken as the propellant depletion 
time, although greater burn times would also produce acceptable re- 
entry. However, for the service module from some abort times, iso- 
lation on the desired re-entry window conditions could not be obtained 
for any vehicle burn time. As discussed previously, when this became 
apparent, the re-entry flight path was changed from 94° to 99*5° • 
Although this new value permitted successful trajectories from all 
abort times, for some abort times the maximum burn time achievable was 
less than the service module maximum burn time. 

An example of the results from this part of the analysis 
is illustrated by Figure 3* ^'or abort from a time of 503 sec. on the 
Boundary Boost Trajectory A, the data show the re-entry velocity and 
associated range angle for burn times from the minimum to the maximum 
for the S-1VB stage. For these data, range angle is unspecified and 
the COV analysis produces a maximum and a minimum velocity for each 
burn time. Both of these velocities are shown as well as the 0 
associated with each. Similar data for Boundary Boost Trajectory B 
are shown in Figure 4. 

3. The final step consisted of incorporating an OC term 

with the analysis as a parameter to be varied to determine the ex- 

treme range angles for a select abort stage burn time. This amounts 
to using the curves of Figures 3 and 4 as a starting point and varying 

c>C o until maximum and minimum ranges are obtained. After carrying 

out the work for several abort times, it became apparent that for a 
given abort time the maximum range was obtained by performing this 
0 incrementation technique with the maximum burn time point (point 
(a) of Figure 3) and that the minimum range was obtained by incre- 
menting the minimum burn time point (point (b) of Figure 3)* Per- 
forming this procedure for a select abort time on boost Trajectories 
A and B provided the maximum and minimum range angles for that abort 
time. 


RESULTS 

As previously stated, only two boost trajectories have been 
used in this abort study. Therefore, for each abort time there exists 
two maximum and two minimum ranges, i.e., one for each trajectory. 
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The results presented below define the greater of the two maxima as 
the maximum and the smaller of the two minima as minimum range. 

For the service module, Figure 5 gives the maximum and 
minimum ranges as a function of abort time. The abort times shown 
represent abort from sub-orbital conditions on the boost trajectory. 
After the curves were calculated, an evaluation of the resulting re- 
entry velocities and flight path angles was made to determine if the 
re-entries were "safe". A safe re-entry is defined to be one which 
does not exceed a pilot acceleration dose limit or one which does not 
skip out of the earth’s atmosphere after re-entry initiation. (See 
Raytheon's "Re-Entry Corridor for Manned Lifting Vehicle" elsewhere 
in this report, for re-entry limits). The curves of Figure 5 were 
completed before these limits were available and, therefore, are 
independent of such constraints. It was found that for the maximum 
range curve, acceleration dose limit was exceeded for abort times 
greater than 606 seconds. For times greater than 606 seconds, the 
maximum permissible 0 would be determined by the acceleration dose 
limit constraint. It is noteworthy that the limits are extremely 
sensitive to re-entry angles. For example, changing re-entry angle 
to 97° would result in exceeding neither an acceleration limit nor 
a skip-out limit. 

Similar information for the sub-orbital abort with the 
S-IV B stage is given by Figure 6. The advantage of the higher thrust 
to weight of the S-IV B compared to the Service Module is evidenced by 
the greater re-entry range capability. As in Figure 5, the curves of 
Figure 6 are independent of acceleration or skip-out constraints. 

Later evaluation showed that for abort times greater than 352 seconds 
the maximum range would be determined by the skip-out constraint. 

For service module abort from super-orbital conditions, 
Figure 7 gives the range capability. Super-orbital abort occurs 
after leaving the waiting orbit. However, the range angles shown 
are for powered flight only, i.e., 0 during coast in orbit is assumed 
zero. The intersection of the two curves represents the maximum 
flight time from which no-coast abort with the service module is 
possible. This time is somewhat less than the full Saturn V burn 
time so that during the last few seconds abort cannot be accomplished 
with the service module. This conclusion is modified if coasting 
were permitted, if a different re-entry angle were used, or if a 
greater thrust to weight ratio were available. 

CONCLUSIONS 

1. The data presented here represents the abort no coast 
range capability envelopes of the Service Module and S-IVB stage from 
the selected boost envelope "boundary" trajectories with specified 
re-entry conditions. 
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2. Some desired re-entry conditions could not be reached 
when the SM was used as the abort vehicle for some booster flight 
times due to its low thrust to weight ratio. 

3. Abort with the SM is impossible after a flight time of 
938 seconds under the ground rules of this study. 

4 

4. It may be possible to extend the 2-dimensional range 
capability by investigating the other values of re-entry angle . 

Other re-entry angles also have a marked effect on range limits de- 
fined by excessive acceleration dose and skip-out. 

RECOMMENDATIONS 

The information of this report should be used to make a 
preliminary selection of desirable abort site locations along the 
AMR boost ground track. Once such sites are chosen, volumes of 
abort trajectories may then be generated from which the powered 
abort guidance equations may be produced. 

Abort using coast and abort into a parking orbit should be 
considered as methods of extending the limiting abort time of the 
SM as established by this study. 
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TABLE 1 

Stage I 

Sea Level Thrust (lbs.) 

Sea Level Specific Impulse (sec) 

Lift off Weight (lbs) 

Propellants consumed (lbs) 

Weight dropped at S-1C separation (lbs) 
Azimuth Angle at Lift-off (deg.) 

Stage II 

Vacuum Thrust (lbs) 

Vacuum Specific Impulse (sec) 

Lift-off weight (lbs) 

Propellants consumed (lbs) 

Weight dropped at S-II separation (lbs) 

Stage III 
Vacuum Thrust 

Vacuum Specific Impulse (sec) 

Lift-off Weight (lbs) 

Propellants consumed to orbit (lbs) 
Weight lost in orbit (lbs) 

Propellants consumed to injection (lbs) 
Total propellants consumed (lbs) 

Weight dropped at S-IVB separation (lbs) 

Gross Payload (lbs) 

* Specific impulse is classified. 


7,500,000 

* 

6,000,000 

4,242,362 

391,560 

70 


1,000,000 

* 

1,366,078 

919,010 

87,800 


200,000 

* 

359,268 

78,521 

5,000 

146,831 

225,352 

23,100 

105,815 
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Figure 2 



I SATURN Y. POWERED ABORT 1 

S-JE B 


ABORT FROM UPPER S -3E B 

BOUNDARY TRAJECTORY A F * 100,000 LBS. 
ABORT TIME ~ 503 SECS. W 0 = 35 9,£68 LBS. 



Figure 3 


84 


| SATURN 31 POWERED ABORT 

S-IZB 


ABORT FROM LOWER 5 - 12! B 

BOUNDARY TRAJECTORY B F= 200,000 LBS . 
ABORT TIME = 503 SECS. W 0 - 359,263 LBS. 



Figure 4 
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SATURN! g POWERED ABORT 

SERVICE MODULE 


SERVICE MODULE 
F = £1,500 LBS. 

W„ = 61, I &0 LBS. 


RE-ENTRY FLIGHT 
PATH ANGLE = 99.5 DEGS. 
RE-ENTRY 
ALTITUDE = lao KM. 



Figure 5 
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SATURN ZZ POWERED ABORT 


S - 12! B 


S-IEB 

F * ZOO ,000 LBS. 
W 0 =: 35 9, 268 LBS. 


RE-ENTRY FLIGHT 
PATH ANGLE = 94 0 
RE- ENTRY 
ALTITUDE = 120 KM . 



Figure 6 


RE-ENTRY RANGE ANGLE ^ DEGS. 
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SATURN g POWERED ABORT 

SERVICE MODULE 


SERVICE MODULE RE-ENTRY FLIGHT 

F = £1,500 LBS- PATH ANGLE * 94° 

W c = Gi,iao LBS- 



Figure 7 
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^ Summary 

n £* I. * * * 5 A 

A study has been made to determine the extreme re-entry 
flight-path angles, for various re-entry speeds, which permit safe, 
unpowered descent for a given manned lifting vehicle. Circular arcs, 
called entry arcs, were located at the initial altitude of 120 km. such 
that, for a specified initial speed and flight- path angle, entry at any 
point within the associated arc ensures the ability to arrive at the des- 
ignated target on the surface of the earth. . - 


I. INTRODUCTION 


This study is part of the general abort re-entry problem. The 

trajectory of a s # pace vehicle may include a boost phase which trans- 
ports the vehicle through the atmosphere and into space, a powered 

and/or cruise phase through space, and a re-entry into the atmosphere 

followed by an unpowered descent to an altitude at which the landing 
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phase can begin. Somewhere along the powered phase of the trajectory 
a decision to abort may be necessary. The capability of a manned ve- 
hicle after an abort, especially its ability to re-enter the atmosphere 
safely, is of utmost importance. 

This report is concerned with the re-entry corridor of a 
manned lifting vehicle descending through the atmosphere. Definition 
of this corridor requires finding the bounds on the re-entry velocity 
and position which permit a safe descent to a specified landing site. 

The acceptable initial conditions for re-entry and descent prescribe 
the terminal conditions for the exo- atmospheric phase of the trajectory. 
The results of the studies of the spatial and atmospheric phases of the 
trajectory must be combined to determine the situations from which 
successful aborts may be initiated. 

II. LIST OF SYMBOLS 


C D drag coefficient 

C lift coefficient 

L 

D drag 

L lift 

P pilot penalty function 

R earth radius 

S reference area 

T final time 

V vehicle speed 

initial (re-entry) vehicle speed 

a aerodynamic acceleration 

g gravitational acceleration 

g g at surface of the earth 

o 

h altitude 

m vehicle mass 

t time 

a control variable, angle of attack 

0 flight-path angle 

0 Q initial (re-entry) flight- path angle 

p density of air 

T pilot acceleration- endurance time 

0 angular displacement, range 
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III. PROBLEM DESCRIPTION 


A. Dynamic Model 

The vehicle is considered to be a particle of constant mass 
which moves in a plane with respect to a spherical, non- rotating 
earth. It is subject to the action of three forces: the inverse- square 
gravitational field of the earth, its lift, and its drag. The descent of 
the vehicle is controlled by varying the lift and drag forces. The 
variation of air density with respect to altitude is included. The force 
diagram is shown in Figure 1. 

The weight of the vehicle is 8500 pounds and its reference area 
is 12. 97 square meters. The aerodynamic coefficients are functions 
of the control variable, a, as shown in Figure 2. The maximum lift- 
to-drag ratio is 0.82 which occurs at a = 50°. During this study, the 
angle of attack was constrained to the interval of -70° to 70° because 
this interval includes the extreme variations in lift, drag, and lift-to- 
drag ratio. 

Below Mach 2, the aerodynamic coefficients are functions not 
only of a but also of Mach number. In view of this consideration, a 
speed of Mach 2 served as the stopping condition for the computing. 
This speed occurs at altitudes compatible with the initiation of the 
landing phase. 


B. Equations of Motion 


where 


mV = - D - mg cos 0 
mV (e + (j>) = - L + mg sin 0 
(R + h) (J> = V sin 0 
h = V cos 0 


g = g. 


R 

R + h 


D = j C D (a) p(h) V 2 S 


L = \ C L ( tt ) p(h) V 2 S 
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and 


p(h) is given by ARDC Model Atmosphere 1956. 
Cj_ ) (q') and C^(a') are shown in Figure 2. 


/ 2 

g Q = 9. 815 m./sec 
R = 6. 371 x 10 6 m . 
S = 12. 97 m. Z 



C. Pilot Acceleration- Endurance Constraint 


For a manned re-entry, a programs which produce excessive 


aerodynamic accelerations must be excluded. This condition is im- 
posed during the solution procedure in the following way. A man's 
ability to remain usefully conscious is a function of both the aero- 
dynamic accelerations he experiences and their durations. It has 
been shown that he can tolerate quite high accelerations if they are 
sufficiently brief. The dimensionless aerodynamic acceleration, a, 
is defined by 


Experiments have yielded the endurance limit r(a) of experienced test 
pilots to given aerodynamic accelerations. By adding the equation 


to the equations of motion, the "acceleration dose" or terminal value 
of the "pilot penalty function" is given by 


where T is the time of flight. When this quantity becomes 1, the 
pilot is assumed to have had a full dose of acceleration; therefore, 
he should not be exposed to further accelerations that would increase 
this dose, if he is to function usefully. Thus, a terminal constraint is 



T(a) 


T 



0 


P<; 1. 
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The pilot acceleration- endurance function, r(a), used in this 
study is shown in Figure 3. This function was derived principally 
from information in References 1-3. The more recent data in 
References 4 and 5 reveal that the function of Figure 3 is conserva- 
tive by factors from 2 to 5, in terms of permissible time for a given 
acceleration, if the pilot is oriented in the most favorable attitude. 

In the current study, however, the attitude of the vehicle is subject to 
wide variations in some maneuvers. If the pilot is exposed to similar 
variations in attitude, he may experience situations where, according 
to Reference 5, his endurance is significantly less than that shown in 
Figure 3. It is believed, however, that the T(a) relation employed in 
the current study represents a reasonable compromise for the speci- 
fication of pilot endurance to acceleration. As a refinement in the 
future, acceleration endurance might be introduced as a function of 
both aerodynamic acceleration and pilot attitude. 


D. Procedure for Trajectory Optimization 


The differential equations of motion together with the pilot 
penalty function form a non-linear system. The initial conditions are 
the altitude, 120 km., the re-entry speed, and the flight-path angle, 
the latter two being parameters of the study. The terminal condition 
is that the pilot not receive more than a full !, dose n of acceleration, 
i. e. P(T)^1. The control variable, a(t), occurs as an unspecified 
function. 

For a given re-entry velocity (speed and flight path angle), 
the end points of the entry arcs are found from the maximum and 
minimum range trajectories. Several nominal a programs are as- 
sumed and the system is integrated by means of a high-speed digital 
computer for each one in turn. None will, in general, yield the ex- 
tremal range nor satisfy the penalty constraint. The most promising 
a program is chosen and then subjected to successive improvements. 
A systematic procedure for producing such changes is the steepest- 
ascent method developed at Raytheon Company and described in 
Reference 6. It is a calculus of variations technique and alters the a 
program in such a way as both to improve the n pay-off !f quantity 
that is being extremalized and to meet any terminal constraints. 

Thus a sequence of a programs is generated. This procedure is 
terminated when the terminal constraints are satisfied, and negligible 
gains in the pay-off quantity are produced by successive iterations. 

For a given re-entry speed, the maximum re-entry flight- path 
angle can be found by increasing the initial flight-path angle until it 
becomes impossible to constrain the pilot penalty function to 1 either 
as a terminal constraint in an extremal- range series of iterations or, 
in more difficult cases, as the pay-off quantity in a minimal- pilot- 
penalty series of iterations. The minimum re-entry flight- path angle 


J 
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is established theoretically at 90° for subcircular entry speeds and 
at an angle that exceeds 90° by an arbitrarily small amount for cir- 
cular entry speed. For a supercircular entry speed, the minimum 
re-entry angle can be estimated using the work of References 7 and 
8, and verified by decreasing the initial flight- path angle until, using 
maximum negative lift, the vehicle rises above the specified maxi- 
mum altitude following the initial pass through the atmosphere. 


IV. RESULTS 

The results are tabulated in Table 1. For a given entry velo- 
city, the angular distance between the target and the re-entry point 
nearest to the target is indicated in the Minimum Range column. If 
the vehicle enters the atmosphere at this distance from the target, 
the descent must be made using the a program associated with the 
minimum- range trajectory. An entry closer to the target will cause 
overshoot because the steepness of the trajectory is limited by the 
pilot penalty function. Similarly, the numbers in the Maximum Range 
column indicate the farthest from the target that entry may occur. 

The entry arc is the circular arc at the specified initial altitude of 
120 km. joining the nearest and farthest possible re-entry points. 
Entry at any point within this arc with the associated initial speed and 
flight- path angle ensures the ability to arrive at the target. The 
entry flight- path angle, which is the direction of the initial velocity 
vector measured counterclockwise from the local vertical, can be 
confined to lie between 90° and 180°. The trajectory for an entry 
flight-path angle lying between 180° and 270° is the same as for its 
mirror image in the 90° to 180° range. 

Table 1. Tabulation of Results 


Entry- 

Entry Flight- 

Minimum 

Maximum 

Entry- Arc 

Speed 

Path Angle 

Range 

Range 

Length 

(m./sec. ) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

750 

90 

0. 93 

1. 05 

0. 12 


180 

0 

±0. 05 

0. 10 

3500 

90 

4. 8 

7. 3 

2. 5 


no 

2. 3 

3. 1 

0. 8 

7833 

90. 50 

46 

138 

92 


101. 75 

8 

16 

8 

11080 

94. 71 

22 

. ❖ 
160, 

138 


99. 8 

15 

no* 
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❖See discussion in text. 
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The re-entry corridor, as it appears in the initial- flight- path- 
angle, initial- speed plane, is shown graphically in Figure 4. 

For the lowest entry speed studied, 750 m./sec., there is no 
restriction on the initial flight-path angle. The vehicle can enter the 
atmosphere with a horizontal velocity or one which is straight down, 
but the range and entry- arc length are so small as to be negligible 
when compared with the performance at higher speeds. 

The steepest entry angle for an entry speed of 3500 m./sec. 
is approximately 110°. At this entry angle, the pilot penalty con- 
straint can be held to 1 for maximum and minimum ranges through 
appropriate modulation of the a program. A critical search was not 
made to verify the possibility of steeper entries because available 
information concerning the entire abort-trajectory problem indicated 
that re-entries for initial speeds of roughly 3500 m./sec. most likely 
will occur for angles less than 110°. The range and entry- arc capa- 
bilities at this speed may be of some significance for an entry at an 
angle of 90° but they both decrease drastically as the entry angle be- 
comes steeper. 

For true circular entry speed, the shallowest possible entry 
angle is undefined. A horizontal circular velocity, 0 o =9O°, results 
in a circular orbit and consequently no entry if the effects of aero- 
dynamic drag are absent. Any initial flight-path angle greater than 
90° will result in re-entry, and the closer this angle is to 90°, the 
larger the maximum range. Similarly, a slight reduction in initial 
speed and/or the presence of slight aerodynamic drag at the specified 
initial altitude will lead to entry. For the solutions obtained during 
this study, the initial speed was circular for the entry altitude, but 
the atmospheric density, and hence drag, were defined to above this 
altitude in accordance with the ARDC Model Atmosphere, 1956. 

In the circular- speed- entry studies, an arbitrarily selected 
shallow initial angle of 90. 50° was found to lead to a maximum range 
of only 138°. As the entry angle becomes steeper, the maximum 
range decreases until it is only 16° for steepest permissible entry, 
6 0 = 101. 75°. The minimum ranges and the entry- arc lengths also 
are markedly less for the steeper entry angles. This situation is 
illustrated in Figure 5. It is significant to note that the entry arcs 
for the extreme entry angles do not overlap; consequently, several 
target areas will be necessary to effect successful recovery of space 
vehicles re-entering at circular speed if initial flight-path angles lie 
anywhere between the limits of 90° and 101. 75°. 

When a vehicle travelling at supercircular speed re-enters 
the atmosphere at a shallow flight- path angle, the aerodynamically 
produced deceleration may be insufficient to prevent the vehicle from 
rising above a specified altitude limit. Thus, the re-entry problem 
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reduces to the determination of the shallowest initial angle that leads 
to the satisfaction of the altitude restriction when the vehicle is flown 
with maximum negative lift. Through the use of results of theoretical 
analyses, as verified by numerical solutions, it was established that, 
for an entry speed of 11, 080 m./sec. (essentially escape speed), ac- 
ceptable re-entry can be accomplished for an entry angle as shallow 
as 94. 71°, but not for one of 94. 55°, when the altitude limit is 
150 km. In lieu of attempting to define 0 Q more exactly within this 
narrow range, 94. 71° was taken as the shallowest initial flight-path 
angle at this speed. 

The steepest entry angle at escape speed is limited by the 
pilot acceleration dose during the initial dive into the atmosphere. 

This dose is critically dependent on the precise modulation of the 
angle- of- attack program. For an entry angle of 99. 8° an acceptable 
pilot- penalty value was achieved for both minimum and maximum 
range trajectories. Among the many trajectories evolved during the 
study of performance for steeper initial flight-path angles, none 
yielded an acceleration dose as low as 1. 

Minimum- range capability for escape- velocity entries also 
is limited by the pilot acceleration dose. For the entry angles 
studied, this range decreased from 22° for 94.71° to 15° for 99.8° 

In the case of the shallow entry angle of 94. 71°, a sufficient margin 
of negative lift was available to prevent the minimum- range trajectory 
from leaving the atmosphere following initial entry. Of course, for 
the actual shallowest permissible entry angle, which is between 94. 55° 
and 94. 71°, the minimum- range trajectory would include a rise to 
the specified maximum altitude of 150 km and the resulting range 
would be substantially greater than 22°. 

The computation of the maximum range for entries at escape 
speed becomes particularly difficult as the steepness of the entry 
angle increases. In these situations, the angle- of- attack program 
during the first 10% or less of the total flight time must be modu- 
lated extremely accurately in such a way that both the pilot- penalty 
and maximum- altitude restrictions are satisfied in a manner com- 
patible with maximization of the range. The total pilot penalty is 
realized during roughly 2% of the flight time shortly after initial 
entry into the atmosphere, and the maximum altitude restriction, 

150 km., occurs later in the flight during a long interval when the 
aerodynamic forces are negligible, thus complicating the solution 
process. During this study, the range capability was computed both 
by optimizing the performance during the entire time of flight and by 
combining extremal solutions for appropriately defined portions of 
the over-all trajectory. Cross checks were made to establish the 
compatibility of these approaches and to ensure the relative validity 
of the answers. The maximum ranges given in Table 1 represent 
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the "best" answers obtained. These ranges definitely are realizable 
under the specified conditions and perhaps can be increased through 
appropriate changes in the angle- of- attack program early in the 
flight. 


The entry arcs for escape- speed entry are shown in Figure 6. 
For entry angles of 94. 71° and 99. 8°, the entry arcs overlap to a 
large extent indicating the feasibility of using a single recovery area. 

In aborts during space missions, of course, the re-entry 
velocities are not subject to close control; they will lie between broad 
limits which are determined by many factors. Based on the results 
given in Table 1, if the speeds may be anywhere in range from zero 
up to escape and entry flight- path angles are unrestricted, recovery 
facilities would have to be provided on a continuous basis throughout 
possible re-entry areas. As the range of expected speeds decreases, 
and as probable flight-path angles are defined, projections may be 
made as to the discrete number of landing sites needed to effect suc- 
cessful recovery. 


V. CONCLUSIONS 


As re-entry speeds increase from 750 m./sec to escape 
speeds, restrictions arise on the possible re-entry flight-path angles. 
The shallowness of the entry, for supercircular entry speeds, is 
limited by the tendency of the vehicle to skip out; the steepness of 
the entry for all except the lowest speeds, by the acceleration- dose 
constraint. Stringent restrictions on initial flight-path angles which 
occur for escape- speed entries are coupled with wide tolerances on 
re-entry position. Re-entry speeds and flight-path angles must be 
limited more than indicated by the results reported here if a small 
number of landing sites is to offer a high probability of successfully 
recovering aborted spacecraft. 

The results of this study define extreme re-entry conditions 
for the specified vehicle when subject only to pilot- acceleration- dose 
and altitude constraints. The re-entry corridor may be changed if 
any of the following considerations are included: the total heat and/or 
heating rate is constrained; the pilot- acceleration- endurance function 
includes pilot attitude dependence; or the magnitude of the angle of 
attack is limited. 
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